Nonclassicality of pure two-qutrit entangled states 
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We report an exhaustive numerical analysis of violations of local realism by two qutrits in all 
possible pure entangled states. In Bell type experiments we allow any pairs of local unitary U(3) 
transformations to define the measurement bases. Surprisingly, Schmidt rank-2 states, resembling 
pairs of maximally entangled qubits, lead to the most noise-robust violations of local realism. The 
phenomenon seems to be even more pronounced for four and five dimensional systems, for which we 
tested a few interesting examples. 
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I. INTRODUCTION 

John Bell [l| has shown that local realistic models are 
impossible for quantum mechanics of two qubits. Af- 
ter some years researchers started to ask questions about 
the Bell theorem for more complicated systems (com- 
pare the classic review in which such problems are not 
covered, Q). First results concerning possible exten- 
sions to entangled states of pairs of d-state systems, with 
d > 3, arrived in early eighties f3[. The early research 
was confined to Stern-Gerlach type measurements (which 
are parametrized by only two numbers, defining the di- 
rection of the measurement axis). Some authors spec- 
ulated that correspondence principle suggests that non- 
classicality should diminish with growing d. However, on 
one hand in the meantime we witnessed a surprise in the 
form of the GHZ theorem for three or more qubits 
the conflict between local realism and quantum mechan- 
ics is much sharper than for two. Thus, with the increase 
of the dimensionality of the full system, non-classicality 
of quantum correlations can grow. On the other hand, 
Peres and Gisin [s'] considered dichotomic observables ap- 
plied to maximally entangled pairs of qudits, and showed 
that violations of the CHSH inequalities robustly survive 
in the limit of d — >■ cx). 

The concept of multiport interferometers, first dis- 
cussed in the context of quantum entanglement by 
Klyshko fe^, gave a hope for operational realizations of 
unitary transformations much richer than those linking 
measurement Stern-Gerlach bases for higher spins. Pro- 
posals of Bell experiments with the multiports were pre- 
sented in 0, II]. Multiport interferometers were shown 
to be capable to reproduce all finite dimensional uni- 
tary transformations 0, i.e., one can have access to the 
fidl U(d) group. This lead to the discovery that two 
maximally entangled qudits violate local realism more 
strongly than qubits, and that this violation grows with 
d, see The strength of the non-classicality was 

measured via a "white" noise resistance. One takes 
a family of mixed states g = VQstate + (1 - v)Qnoise, 
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where Qstate represents the state under consideration, and 
Qnoise = ^ 1 IS the maximally mixed state. The lower 
is the threshold value v = Vcrit beyond which state does 
not violate local realism in a given Bell experiment, the 
bigger is the noise admixture 1 — Vcrit , which is needed to 
erase the non-classicality of the state in the experiment. 
Note that Vcnt is a natural generalization of the "factor 
of violation" of a Bell inequality. E.g. for the GHSH in- 
equality, if B is the Bell operator, its local realistic bound 
can be violated maximally by a factor oivZ}^^ . The same 
can be said about CGLMP inequalities [Tl|. The results 
were obtained via a numerical analysis, which is a proto- 
type of the one which will be presented. The values were 
confirmed by Collins et al. [11| who derived a set of tight 
correlation Bell (CGLMP) inequalities specific for qudit 
measurements. This tool helped to discover a strange 
property of two qudit states. Maximally entangled states 
do not violate the CGLMP inequalities for qutrits maxi- 
mally. The Schmidt decompositions of the optimal states 
does have all amplitudes of the same modulus — one of 
them is smaller [l3|. The results were generalized fur- 
ther by Chen at al. With a new version of a qudit 
Bell inequality, the optimality of non-maximally entan- 
gled states was shown in arbitrary dimension, |15|. Note 
however, that in the papers a specific set of local observ- 
ables was used, defined by a set of local phase shifts (in 
the Schmidt decomposition basis for the state), followed 
by an "unbiased" multiport beamsplitters (class Ml in 
Fig. 1), and finally detectors. Thus, non-classicality was 
not fully mapped. 

Here we report an exhaustive numerical analysis of the 
two-qutrit Bell experiment involving two settings per ob- 
server. The observables are absolutely general, all possi- 
ble ones are studied (that is the unitary transformations 
defining the measurement bases are forming the U(3) 
group). Such an analysis allows us to map the strength 
of violations of local realism (as measured by noise re- 
sistance) for all possible pure two-qutrit states. In this 
way we discover that Schmidt rank-2 states give rise to 
highest violations. This phenomenon continues if we in- 
crease the number of settings for the observers to three 
and four. The investigations were also extended to higher 
dimensional systems, for which we tested a few specially 
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chosen states. The phenomenon seems to persist, and is 
even more pronounced. As a by-product of our investiga- 
tions we have found for two qutrits a relatively simple set 
of unitary transformations (defined by just three phases), 
which gives a map of violations of local realism which di- 
verges from the one for the full U(3) group defined mea- 
surements by maximally just 1.5%, and for a broad class 
of states leads to the same level of non-classicality. We 
have also studied some interesting mixed states, which 
contain bound entanglement. 

II. DESCRIPTION OF THE METHOD 

Let us move on to the quantum processes that we 
analyze, and the numerical approach. We begin with 
the two-qutrit case. In our numerical analysis (called 
stesun-roller) we consider a class of pure states of two 
qutrits, which contains all possible Schmidt decomposi- 
tions: 

\ip{a,(3)) = cosa|00) -|-sina(cos/3|ll) -f sin;3|22)). (1) 

Two spatially separated observers perform measurements 
of m alternative local noncommuting trichotomic observ- 
ables: Ai,A2, ■■■Am for Alice and Bi,B2, ...-Bm for Bob. 
We assume that they measure observables defined by a 
set of phase-shifters and (see Fig. [1]): 

• (Ml) one unbiased three-port beam-splitter (frit- 
ter, for its properties, see|8[), 

• (M2) two unbiased three-port beam-splitters, 

• (M3) three unbiased three-port beam-splitters, 
and in the ultimate case 

• (U(3)) any three-dimensional unitary transforma- 
tion belongs to the U(3) group (for parametrization 
see e.g. [l3|)- 
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FIG. 1. Measurement devices represented by unitary trans- 
formations Ml, M2 and M3. The inputs are on the right hand 
side. Boxes represent unbiased symmetric three-port beam- 
splitters 

A d input and output unbiased multiport performs a 
unitary transformation described by the Fourier matrix 
Uki = rf-i/2e«(2fei7r/d)^ kj = Q^ l...,d- 1. The uni- 

tary transformations Ml and M2, obviously cannot re- 
produce the full U(3) group. We have checked that this 
is also the case for M3 transformations. 



By saying that an experiment is local realistic we un- 
derstand that it has a local realistic model for the as- 
sumed set of settings. In order to obtain the value of the 
critical visibility Vcrit to allow such models for any ob- 
servables from a given class, we follow the procedure of 
the kind first used in Refs. [13, [TB] (but with a much more 
advanced implementation and structure). The method is 
described in the Appendix. 



III. RESULTS 

We have computed a map of strength of violation of 
local realism by pure two qutrit states. There one can 
find a range of the a and /3 parameters of the state 
([1]), for which the critical visibility is below the low- 
est known critical visibility (0.6861) in the case of the 
CGLMP inequality, for the asymmetric state \ipasym) 
found in ref. [13] (see Fig. [2]). The lowest critical visibil- 
ity is 0.6821 and surprisingly corresponds to the Schmidt 
rank-2 states Wsym'"^)- |V'(90°, 45°)) = (|11)-H|22))/V2, 
or |V'(45°,0°)) = (|00) + |11))/V2 and |i/'(45°, 90°)) = 
(|00) + |22))/^/2 (outside the map, but trivially related 
with the previous one) . We can observe this phenomenon 
by using the most general local three-dimensional observ- 
ables (U(3)) and transformations M3, while for less gen- 
eral observables (e.g. Ml, M2), it is unobservable. 

Despite the fact that the model noise well describes 
e.g. dark counts, one may criticize "white" noise resis- 
tance as a measure of non-classicality, see e.g. [13] ■ Thus, 
we have studied also other types of noise. For a prod- 
uct noise admixture, of the form pA ® Pb, where pa(b) 
is a reduced density matrix of the state ([T]) of the sys- 
tem A(B), the critical visibility for the state IV'Iym'^^^) is 
equal to 0.7071, whereas for the states \4'asym), IV'sym) 
is the same as in the case of white noise admixture, 
respectively 0.6962, 0.6861. For a dephasing noise, 
cos2al00)(00l-t-sin2a(cos2/3111)(lll-|-sin2/3l22)(22l), we 
observe violations of local realism for any state \il]{a,l3)) 
(a > 0) and any amount of the noise {v < 1). Such a 
feature is well known for qubits, however for qutrits thus 
far it was only a conjecture. 

We can also observe such phenomena in the 4- 
dimensional case. For the most general local 4- 
dimensional observables defined by the U(4) group, the 
critical visibility obtained for a symmetric Schmidt rank- 
2 state (e.g. (JOO) + \n))/V^) is equal to 0.6442, whereas 
the critical visibility necessary for violation the CGLMP 
inequality is equal to 0.6728 and 0.6906 for the asymmet- 
ric [14|] and symmetric ll] states respectively (see Tab. 

The critical values for rank-2 states. Tab. HI agree 
with those reported in [T2| . where such numbers were 
found employing an ad hoc approach, aimed at showing 
conceptual problems related with the noise resistance as 
measure of non-classicality. Here we establish that these 
are the lowest possible critical visibilities for two-qutrit 
states and give a strong numerically supported conjee- 
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0.6071 



TABLE I. The critical visibilities for special states. The symbol \'^^sy^'^) for k = d stands for a maximally entangled state of 
two qudits, that is with symmetric (all equal) Schmidt coefficients. If fc < d it is a state of two qudits resembling a maximally 
entangled symmetric state for fc-dimensional systems. The states ItpasyrrT'') represent in the case k = d the one which maximally 
violates the CGLMP inequality for qudits, whereas for k < d, they are two qudit states resembling the one optimal for the 
CGLMP inequality for a k dimensional problem. The observers perform unrestricted measurements determined locally by U(d). 
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tage of the numerical method is even better visible and 
occurs already for a > 70° (see Fig. [^b). 
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FIG. 2. The map of critical visibilities Vcrit, that is noise re- 
sistance of non-classical correlations for the state |-i/)(a,/3)) — 
cosa|00) -I- sin a(cos /3|11) -f- sin/J|22)) for the case in which 
both observers use any observables. li v > Vcru, there does 
not exist any local realistic model describing quantum prob- 
abilities of experimental events. 



ture that it is so for higher-dimensional systems. One 
could claim that a rank-2 state is effectively a two qubit 
state. However, as it is clear from our graphs, in the case 
of qutrits there is a continuous family of (strongly asym- 
metric) rank-3 states, which violate local realism stronger 
(in terms of noise resistance) than the optimal state for 
violation of the CGLMP inequality. 

We compared our computer code, which is equivalent 
to a full set of Bell inequalities for the given problem, 
with the CGLMP inequality. For two qutrit states of the 
form |-0(q;, 45°)), if we restrict observables to Ml type, 
the numerical method gives the same[? ] critical vis- 



ibility as the CGLMP inequality [llj. Note that such 
observables were used in earlier works, refs (lol - [T5j . In 
the case of observables of M2 kind, for a > 73° the crit- 
ical visibility obtained with stecun-roller is lower than 
predicted by CGLMP inequality (see Fig. |3^). For ob- 
servables parametrized by M3 (and full U(3)) the advan- 
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FIG. 3. The dotted lines represent the critical visibility, 
which is necessary to violate the CGLMP inequality by states 
|'!/'(a, 45°)). The solid lines correspond to the critical visibility 
obtained with steam-roller. M2 (left graph) and M3 (right 
graph) transformations are considered here. The points at 
the curves with an undefined tangent indicate that different 
Bell inequalities are maximally violated by states to the left 
and right of the point. 

For the two-qutrit case, there are some ranges of the 
parameter a in 1-0(0,45°)), for which: 

• the observable generated by M2 gives a better criti- 
cal visibility than the ones obtained for the observ- 
able Ml (a < 49° and a > 73°). 

• the observable for M3 gives a better critical visibil- 
ity than the ones obtained for the observable M2 
(a < 54° and a > 70°). 

All of these cases are presented in Fig. |31 The critical 
visibilities obtained for the M3 and U(3) observables are 



4 



the same for the |V'(ci;, 45°)) state. 

-M3 



/? the highest difference 5vY}it 



For other values of 
(%) between critical 
visibilities obtained for these observables is less than 1%. 
If one increases the number of settings, up to four on each 
side, this does not reduce the critical visibility for any of 
the types of observables. 
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FIG. 4. The critical visibilities for states 45°)). Each 

line corresponds to a different type of the unitary transfor- 
mations (Ml, M2, and M3). Up to four settings per observer 
were used, the relations did not change. 



A. Almost perfect two-qutrit Bell interferometer 

Let us consider the M3 interferometer shown in Fig. [T] 
with (/>! = 05 = 06 = 0- The differences between critical 
visibilities obtained for such transformation and U(3) are 
shown in Fig. [5^ for the state |-0(a, /?)) and in Fig. [Sb 
for the state \ip{a, 45°)). They are less than 1.5% for any 
a and /3, and in many cases the same (including specific 
points related to the symmetric \^psym) and asymmtric 
IV'osj/m) state). Thus if one is interested in situations in 
which violations of local realism play an essential role 
this relatively simple device might be optimal. 



B. Question of violation of local realism by bound 
entangled states 

As we have at our disposal a tool which allows to test 
violation of local realism for any local observables, we also 
analyzed some classes of two-particle bound entangled 
states, namely: the three-dimensional Bennett state [18| . 
the three-dimensional Horodecki state ^19'] and its gener- 
alization to the four- and five-dimensional states [20], the 
four-dimensional Pankowski - Horodecki state [21|. We 
established that for U(3), and in the last case U(4) trans- 
formations (that is all possible local measurements with 
two settings per observer) one does not observe violations 
of local realism. 
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FIG. 5. a) Differences between critical visibilities obtained for 
the U(3) and M3 (with (j>i = 4'5 = 4's = 0) transformations 
for the state \ip{a,j3)). Peak points are lower than 1.5%; b) 
Differences between critical visibilities obtained for the U(3) 
and M3 (with <j>i = (jis = <j>6 = 0) transformations for the 
states |V'(a,45°)). 



IV. FINAL REMARKS 

We have charted unknown territories in the map of 
non-classicality of two qutrit states, and made some ex- 
cursions to higher dimensional problems. The obtained 
results were difficult to anticipate. Our method is quite 
flexible and we plan to move on to the yet uncharted 
multi qudit states. Note, that the computer code, af- 
ter special processing of experimental data is able to 
decide whether the data allow a local realistic model 
or not. Such applications will be found in forthcoming 
manuscripts. The general theory of the code itself, and 
particular solutions, will be presented in ref. [2^ . 
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Appendix: Description of the method 

The critical visibility is calculated with the use of the 
simplex algorithm [23] from the GLPK library (2^ . This 
critical visibility is dependent on the angles (parameters 
of the observables) and it can be shown that it is in fact 
a continuous but non-differentiable function, of the an- 
gles [29| . An intuitive explanation of why this is the case 
will be given in the last paragraph. The critical visibil- 
ity function can be minimized with the use of non-linear 
optimization algorithms. Since it is non-differentiable we 
are limited to the ones which do not require the deriva- 
tive of the optimized function. We use the established 
Nelder Mead method [25|. Despite problems raised in 
[26j where it is shown that the method does not converge 
to a stationary point in the general case, in our case 
it converges nicely when applied to the critical visibil- 
ity function. This was tested extensively by comparison 
with both results derived analytically and known from 
other sources. It should be noted, however, that while 
the search for critical visibility is guaranteed to return 
the correct result 23] , the Nelder- Mead method may re- 
turn a local minimum. Therefore, one can be sure that 
the actual minimal critical visibility is not greater than 
the one given here (but it might be smaller). 

For qutrit experiments with two alternative measure- 
ment settings per observer {m — 2), realistic models are 
equivalent to the existence of a joint probability distri- 
bution p;,.(ai, 02, &i, 62)7 where = 0,1,2 (i = 1,2) de- 
notes the result of the measurement of Alice's i-th ob- 
servable (Bob's results are denoted by bk = 0,1,2 {k — 
1,2)). Quantum predictions for the probabilities should 
be given, if the model exists, by marginal sums 

2 

P{ai,bk\Ai,Bk) ^ ^ Pir{ai,a2,bi,b2), (A.l) 

where P{ai, bk\Ai, Bk) denotes the probability of obtain- 
ing the result by Alice and bk by Bob, if they measure 
the observables Ai and Bk, respectively and i' ^ i,k' ^ k. 

If we admix some amount of white noise to the 
two qutrit state, we obtain the quantum probabilities 
PQM{ai,bk\Ai,Bk) = P{ai, bk\Ai, Bk) in the form of 

P{ai,bk\A„Bk) = vPQM{a„bk\A„Bk) + ^{l-v), 

(A.2) 



where Pqm denotes the quantum-mechanically calcu- 
lated probability for the state without the noise admix- 
ture. 

For u = 1, it is known that for some entangled states, 
no local realistic probability distribution pir (ai , 02 , 61 , 62 ) 
which could satisfy the set of equalities (|A.1|) exists. 
Therefore, no local realistic model can reproduce the pre- 
dictions of quantum mechanics for these entangled states. 
However, for such states there always exists the critical 
visibility Vcrit such that for v ^ Vcrit there exists a lo- 
cal realistic probability distribution pir{ai, 02, ^17 ^2) that 
satisfies the set of equalities (jA.ip . Our task — for a 
given state p — is to find settings for which the critical 
visibility is minimal. 

The computation of Vcrit (for the given parameters) is 
equivalent to finding the set of 3** positive numbers (the 
values of pir{cLi, 0-2, bi, 62)) summing up to 1.0 and fulfill- 
ing the 2^ • 3^ = 36 conditions given by (jA.lj) such that 
V is maximal. Therefore, v and Pir{ai,a2,bi,b2) can be 
treated as variables lying in a (3^ + l)-dimensional real 
space. The set of linear conditions (|A.1|) and the con- 
ditions < V < 1;Q < Pi < \ define a convex set in 
this space. Our task is to find a vertex of the set with 
the largest v coordinate. This is done by means of linear 
programming [137 Of course one can easily formulate 
generalization of the presented approach to more settings 
per observer and systems of higher dimensionality. How- 
ever in such cases the computing time increases, and thus 
we studied only specific interesting examples. 

The 36 probabihties P{ai,bk\Ai, Bk) define a point S 
in a 36-dimensional space, which might be inside, on the 
boundary or outside the Pitovsky-Bell polytope in that 
space [271 ]. The point S in this space is within (i.e. in- 
side or on the boundary of) the Pitovsky-Bell polytope 
if and only if the set of equalities (jA.ip is satisfied by the 
coordinates of S. Let us consider a state \S) and a set 
of observables Ai , Bk for which the corresponding point 
S is outside the Bell polytope (for v = 1). The posi- 
tion of S is dependent on v, with S{v = 0) well inside 
the polytope and S{v = Vcrit) on a facet of the poly- 
tope. Given (jA.2|) we know that S{v) is a linear func- 
tion in the 36-dimensional space. The critical visibility 
function for slightly perturbed angles (parameters of the 
observables) can now be calculated by recalculating all 
P{ai,bk\Ai, Bk) probabilities and connecting S{v = 0) 
and the new S'{v = 1) with a straight line. S{v — 0) 
will not change its position since for visibility equal to 
the observables do not matter. S'{v = 1) will be slightly 
perturbed. Thus the new S'{v = Vcrit) point will also be 
"just next to" the previous critical visibility point and 
hence — the critical visibility function is continuous. It 
is also differentiable as long as we stay on the same facet 
of the Bell polytope. But if we perturb the angles in a 
vicinity of an edge of the Bell polytope, differentiability 
is lost. Hence, as mentioned before, the critical visibility 
is a continuous but non-differentiable function of param- 
eters defining the observables. 
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